Time-dependent Maxwell field operators and field energy density for an atom near a 

conducting wall 
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We consider the time evolution ol the electric and magnetic field operators for a two-level atom, 
interacting with the electromagnetic field, placed near an infinite perfectly conducting wall. We 
solve iteratively the Heisenberg equations for the field operators and obtain the electric and mag- 
netic energy density operators around the atom (valid for any initial state). Then we explicitly 
evaluate them for an initial state with the atom in its bare ground state and the field in the vacuum 
state. We show that the results can be physically interpreted as the superposition of the fields 
propagating directly from the atom and the fields reflected on the wall. Relativistic causality in the 
field propagation is discussed. Finally we apply these results to the calculation of the dynamical 
Casimir-Polder interaction energy in the far zone between two atoms when a boundary condition 
such as a conducting wall is present. Magnetic contributions to the interatomic Casimir-Polder 
interaction in the presence of the wall are also considered. We show that, in the limit of large times, 
the known results of the stationary case are recovered. 

PACS numbers: 12.20.Ds, 42.50.Ct 



I. INTRODUCTION 

The existence of zero-point field fluctuations in quan- 
tum electrodynamics has striking consequences such as 
Lamb shift, spontaneous decay, Casimir and Casimir- 
Polder forces, which give evidence of the quantum nature 
of the electromagnetic field [l], 0]. Moreover, Casimir 
forces are a manifestation of the quantum properties of 
the radiation field and of vacuum fluctuations at the 
macroscopic level. Vacuum fluctuations are modified by 
the presence of boundary conditions as well as by atoms 
or molecules. Thus the presence of macroscopic objects 
can change the radiative properties of the atoms. In 
particular, it is known that long-range intermolecular 
forces change when the atoms are in the proximity of 
a boundary condition such as a conducting plate or a 
cavity [H, Q • This is also important from an experimen- 
tal point of view, because a precise comparison between 
measurements of Casimir and Casimir-Polder forces re- 
quires calculating these forces in situations as close as 
possible to realistic laboratory situations. Also, it has 
been recently suggested that appropriate boundary con- 
ditions may enhance Casimir-Polder forces [!, H, @, @], 
and this can be very important for their direct experi- 
mental measurement (usually they are very tiny forces). 

Calculation of time-dependent electric and magnetic 
field operators in the Heisenberg picture generated by 
an atom, even when specific boundary conditions are 
present, are a powerful and general tool to calculate 
various radiative processes for different initial atomic 
states. This kind of calculation have been done some 
years ago by Power and Thirunamachandran for atoms 
in the free space @; they have also successfully applied 



the results for the evaluation of intermolecular Casimir- 
Polder forces, both for ground-state and excited-state 
atoms 0, EH, Cl- 
in this paper, following the same method of Power and 
Thirunamachandran, we calculate the time-dependent 
electric and magnetic Maxwell field operators in the 
Heisenberg representation generated by an atom placed 
near an infinite and perfectly reflecting plate. We will 
also apply our results to the evaluation of the dynami- 
cal (time-dependent) Casimir-Polder interaction between 
two neutral atoms near a conducting wall. We shall 
consider both electric and magnetic interactions. Dy- 
namical Casimir-Polder interactions have been recently 
considered for atoms in the free space also for partially 
dressed states [HI or excited states as well as in 
the case of one single atom and an infinite conducting 
wall [H, [H, [ijj • The method used is very suitable for 
generalization to other initial states of the system. 

This paper is organized as follows. In Section [II] we 
introduce the multipolar coupling model for a two-level 
atom interacting with the radiation field in the electric 
dipole approximation and in the presence of a conduct- 
ing wall. Then we solve the Heisenberg equations for the 
photon creation and annihilation operators up to the sec- 
ond order in the electric charge, using an iterative tech- 
nique. In Sections (TTTJ we evaluate first- and second-order 
corrections to the electric and magnetic field operators, 
discussing their physical properties. Finally, in Sec. IIVI 
we use the results obtained to evaluate the electric and 
magnetic energy density operators and their expectation 
values on the bare atom-field ground state. In Section 
fVl we use the expressions obtained for discussing dy- 
namical electric and magnetic Casimir-Polder interaction 
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energies between two atoms or polarizable bodies when 
a boundary conditions such as the conducting plate is 
present. We remark that the expressions obtained for 
the perturbative expansion of the fields are general and 
do not depend on the specific initial state of the atom- 
field systems. Hence they can be also used for evaluating 
Casimir-Polder energies for different initial states such as 
partially dressed or excited states. 



turn electromagnetic radiation field. We choose a refer- 
ence frame such that the atom is placed in = (0, 0, d) 
whereas the mirror coincides with the plane z — 0. We 
work in the multipolar coupling scheme and within the 
electric dipole approximation [181 [19( . Thus the Hamil- 
tonian describing our system reads 



II. THE HAMILTONIAN MODEL 



We consider a two-level atom in front of an infinite 
and perfectly conducting wall, interacting with the quan- 
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In this expression the field is described by the set of 
bosonic annihilation and creation operators a^j and ■ , 
associated with a photon of frequency cjk = ck, while the 
matrix element of the electric dipole moment operator \i 
and the pseudospin operators S+, S- and S z are asso- 
ciated to the atom, which has a transition frequency luq 



Q. Also, fkj(r) are the field mode functions in the pres- 
ence of the wall, which in Eq. {TJ are evaluated at the 
atomic position ta- Their expressions can be obtained 
from the mode functions of a perfectly conducting cubical 
cavity of volume V = L 3 with walls (-L/2 < x,y < L/2, 
0<z<L) 0,13 



(fkj(r))x = V^ekj)* co 
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(fkj(r)) z = Vs(e kj ) z sin k x (x+ — ^ sin k y (y+—^ cos(k z z) 



where k x = Itt/L, k y = rrnr/L, k z = tvk/L (l,m,n = 
0, 1, 2, . . . ) and ekj are polarization unit vectors. In order 
to switch from the cavity to the wall at z — one has to 
take, at the end of the calculations, the limit L — > oo. 



As discussed in the Introduction, many properties of 
the quantum electromagnetic field and of the atom, in 
the presence of the wall, can be studied in terms of the 
expressions of the time-dependent electric and magnetic 
field operators for the interacting system. In the Heisen- 



berg picture they can be written as 



fkj(r)a k j(*) + h.c. 



V 



kj 



(3) 



V 



where the time dependence is fully included in the pho- 
ton creation and annihilation operators. Their expres- 
sions are formally given by the solutions of the relative 
Heisenberg equations associated to our Hamiltonian op- 
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erator ([I]). However, it is not possible to solve exactly 
these equations for our model. Thus we use an iterative 
technique and the solution for the annihilation operator 
can be written as a power series in the coupling constant 

a kj (t) = aW(t)+a$(t)+aW(t) + ... (4) 

where the contribution a^j (t) is proportional to the i-th 
power of the electric charge. By taking the Hermitian 



conjugate we get the analogous expansion for the cre- 
ation operator. For our purposes we need their expres- 
sions up to the the second order only. The result, which 
is a straightforward generalization to the case with the 
boundary condition of that obtained for an atom in the 
free space [H, HH , is 
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(5) 



where we have introduced the auxiliary function 



III. PERTURB ATIVE EXPANSION OF THE 
FIELD OPERATORS 



d£ 



(6) 



All the operators appearing in the RHS of Eq. (|5|) 
are evaluated at t = 0: from now onwards, when atomic 
and photonic operators are written without an explicit 
time dependence, it is meant they are at t = 0. The 
zeroth-order term does not depend on the presence of 
the source, being only the free-field contribution. As in 
the case of an atom in the free space discussed in [9j, 
the first order term contains only atomic operators, and 
the second order term contains both photonic and atomic 
operators. Substituting the expansions ([4]) and (O into 
©, we get the iterative expansion of the electric and 
magnetic Maxwell field operators 



E(r, t) = E(°' (r, t) + E« (r, *) + E^ (r, t) + 
H(r, t) = H<°> (r, *) + (r, *) + (r, *) 



(7) 



In the next Section we explicitly evaluate the first and 
second order corrections of the fields. 



Using the first of (O , the expression of the zeroth-order 
term is 



E(°)(r,i) 



i\l ^fkj(r)e- Mkt a ki + h.c. 



V 



H<°>(r,t) = ^[V x f k ,(r)]e-^a k , + h.c. 

(8) 



V 



kj 



In order to obtain the first-order correction, we substi- 
tute the second of ([5]) into ([3]), obtaining 

E( 1 )(r,t) = ^^fcf kj (r)( / x-f kj (r^)) 



x (s+F(u + uj k , t) + S_F(w k - w , t)\ e-^ kt + h. 



H (1) (r.«) = ^Et V x f kj -(r)](M-fkiM) 



r 



kj 



x (s + F(u + Wk, t) + S*_F(cj k - w , t)^ e- lWk * + h.c. 

(9) 

We can now take the continuum limit, by allowing 
L — > oo. The continuum limit can be performed with the 
following prescription J2 k j — > rg^p / ^ 3 k- The mode 
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functions f k j (r) are of course affected by this transforma- component of the fields is 
tion since they contain L. The result for the p — x,y, z 



4 13 ( r >*) = X^X>« / rffcfc 3 (5 + F( Wo + Wk ,t) + 5_F( Wk ~ Wo ,i)) e - i ^*^ f rfn(f ki (r)) p (f kj (r A )) 9 
71 q J 3 

H^(r,t) = ^^2 N J dkk^S+F^+u^ + S^F^-u^y-^Y, j «(V x f kj (r)) p (f kj ( rj4 )) f 



(10) 



I 

where the integrations over k run from — oo to +oo, fi q calculated by lengthy straightforward calculations (see 
is the q = x,y, z component of fi. The polarization sum the Appendix for more details). Finally, we get 
and the integral over angular variables can be explicitly 



^ r , n i 



q v=0,l 



q v=0,l 



(11) 



where F pq = 



-V S pq + VpVg, G pq — — e pgs V s (e Q /3 7 be- 



ing the complete antisymmetric tensor), and V p = t^- 
(p = 1,2,3) are differential operators acting on the 
variables x,y,z. The index v = 0, 1 labels two pos- 
sible values of the quantity rv, namely ro = |ro| = 
|r — r^J = ^Jx 2 + y 2 + (z - d) 2 and i\ = |n| = | r — r j | = 
tJx 2 + y 2 + (z + d) 2 , where 17 = (0, 0, — d) is the posi- 
tion of the image of the atom reflected on the wall, ro is 
the distance of the observation point r = (x,y,z) from 
the atom in = (0, 0, d), whereas 17 = (0, 0, —d) is its 
distance from the atom's image. 

The final step for the evaluation of the first-order cor- 
rection is the evaluation of the fc-integral. Using ([6]) and 
the following integral representations of the Heaviside 



step function 0(x) 



1 _ p ±i(k -k)ct 

dk : : sin(fcrv) = Tre ±lk ° r » 6(ct - rv) 



dk 



k — kg 
2 _ pi(fe±fco)ct 



k±h 



i) 



cos(fcrv) = -me ±%k ° r »6{ct - rv) 

(12) 



we obtain the following expressions for the first-order cor- 
rection of the electric and magnetic field operators gen- 
erated by the two-level atom in front of a perfectly con- 
ducting wall 



1 i 

(13) 

ff«(r, t) = -J2 ^hir {1 - S ^^ q G pq -^[S + e tk ^ ct - r ^ + S^e- ik ^ ct -^]9{ct - rv). 

q u=0 Tv 



I 

We wish to stress that expressions (jT3"]l are general opera- elements between any pair of atom-field states. 

torial relations, which can be used for calculating matrix . ...... 

It is worthy analyzing from a physical point of view 
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the expressions obtained. The sum over v in (fT3|) de- 
scribes the existence of two distinct contributions to 
these operators. The first one depends on the parameter 
ro = |r — r^l, that is the distance between the observation 
point r and the atom. We shall refer to it as the direct 
term. The second one is a function of v\ = |r — rj|, 
that is the distance between the observation point and 
the image atom. We will call it the reflected term. Both 
terms propagate satisfying rigorously the causality prin- 
ciple as shown by the presence of the 8 functions. We 
can check that our results (fT5|) correctly yield the known 
analogous results for an atom in the free space as ob- 
tained in Q , if we take the limiting case of the wall very 
distant from the atom: in fact, in this limit ro is finite, 
whereas n grows with the atom-wall distance. In the 
same limit we can thus keep only the v = term in (| 1 3|) , 
and obtain the same expressions found in Q. A final 
remark concerning the action of the differential opera- 
tors in p3p is necessary, because in principle they should 
be applied to the Heaviside functions also. This means 



that the Dirac delta function and its derivatives show up, 
yielding a singular impulse propagating as a wave front 
on the light cone, strongly dependent on the model used 
(in particular on the dipole approximation). Similar sin- 
gular behaviour were found in the field emitted during 
the dynamical dressing of an atom [l3|, El. In the 
applications we will discuss in the following Sections we 
will not consider the electromagnetic field on this (singu- 
lar) wave front and so we will not explicitly evaluate it 
(hence the results in the next Section will not be valid 
for t = ro/c or t = r\/c). 

We now consider the second-order correction to the 
electric and magnetic field operators, obtained by substi- 
tution of the third of ([5|) into |3]). Since the calculation 
is quite straightforward, although lengthy, we shall only 
give the final results. In the expressions obtained, two in- 
dependent sums over wavevectors and polarizations are 
present: one of them can be explicitly performed and fi- 
nally yields a 6 function ensuring relativistic causality. 
Thus we get 



kj q i/=0,l 
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_ /-ik(ct-r v ) _ e ife (ct— r u )\ j \ / g — ife(ct— r„) _ g-ifco(ct-r„)\ 
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_ / g — »fc(ct— Ty) _ c ik (ct-r v )\ | \ l-ik(ct-r v ) _ e -ik a (ct-r u )\ 

^0 + w k ' — " 
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_ /gifc(ct— iy) _ c —ik (ct—r„)\ | 1 / c ik(ct-r v ) _ gifco(ct— r„)\ 



LUQ - CJk 

r 



0(ct-r w ). 



The second-order fields in (Tj"4")) . similarly to the first- 
order ones (fl"3|) , show the presence of a direct and a re- 
flected term. 



IV. ELECTRIC AND MAGNETIC ENERGY 
DENSITIES OF THE FIELD 

In this Section we shall use the results of Sect ion HTT1 for 
calculating the electric and magnetic field energy densi- 
ties generated by the atom when the plate is present; they 
are of course related to the field fluctuations generated 
by the atom. In the next Section the results obtained 
will be applied to evaluate the dynamical Casimir-Polder 
potential between two atoms placed near the conducting 
plate (in the so-called far zone). 



We shall first evaluate the square of the electric field 
operator and its average value on the initial state we are 
mostly interested in. The square of the electric field oper- 
ator coincides with the electric energy density operator, 
apart a factor 1/87T, and also with the electric field fluc- 
tuations, assuming a vanishing expectation value of the 
electric field. Using the expansion of the electric field up 
to the second order, we have 



E 2 {r,t) =E("'(r,i)-E(°I(r,() 

+ (E (Q \r,t) -EW(r,i) +EW(r,<) .E<°>(r,t) 
+ (e^ (r, t) ■ E« (r, f) + E<°> (r, f) • E^ (r, t) 
+ E^(r,<) .E(°)(r,t) 



(15) 
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The first term is a free-field term that does not depend 
on the presence of the source; then we have first- and 
second-order terms. From now onwards we consider the 
specific initial atom-field state |0kj, Ia), where the field 
is in its vacuum state and the "source" atom (A) is in 
its bare ground state, and calculate the (time-dependent) 
expectation value of (|15p . This quantity gives the time 
evolution of the electric energy density during the self- 
dressing process of atom A, that in the next Section will 
be used to evaluate the dynamical Casimir-Poldcr inter- 
action energy between two atoms when the conducting 



plate is present. 

The average value of the zeroth-order contribution is 
independent of the presence of the atom and thus does 
not contribute to atom-atom interactions (it gives the 
time-independent but spatially dependent fluctuations 
of the electric field with the plate present and in the 
absence of the atom). Also, the first-order term has a 
zero mean value on our initial state. Thus, the first rele- 
vant non- vanishing contributions are given by the second- 
order terms. We start evaluating the average value of 
E«(r,i) -E«(r,i) 



(0 kj ,U |E«(r,t).E«(r,i)|O kj ,U) 

= E E M 9 M-ir (1 - 5 ^-ir' (1 - VJ ^ 

pqq' vv'— 0,1 



-6{ct 



F , 



-9{cb-r v ,) 



(16) 



where k = u)q/ c. As discussed in Section IIII1 we do 
not to apply the differential operators to the 8 functions, 
because this would only add singular terms on the light 
cone, in which we are not interested at the moment. Us- 
ing the result 

F pq — = k 3 f pq (kry kr (17) 



where 

/«(H = ^^+(^-^)(^ 3 -5^ 3r ) (18) 

we obtain 



<0 k ,,U lE^M) •E( 1 >(r,i)|0 kj -,U) = - E ^oV^'(-l) ,y(1 ^ ) (-l) ,y ' (1 " VJ / P9 (fco^)/ P9 '(-W) 

x e tko< - r "- r ^0{ct - r v )B{ct - r v .). 



(19) 



We can analyze the previous expression in three different space-time regions. When ct < ro, that is for times shorter 
than the time taken by a light signal to travel from the atom in ya to the observation point r, it is zero because both 
Heaviside functions are vanishing. For ro < ct < rj., i.e. after that a light signal can travel from the atom to the 
observation point but before this point has been reached by the signal reflected on the wall, only the contribution 
with v = v' = is different from zero and the result is the same obtained in [l^] for an atom in the free space. For 
ct > n, that is after that even the reflected signal has passed the observation point, all terms contribute, and we can 
recognize a contribution from the atom [v = v' = 0), one from the image atom (v = v' = 1) and also interference 
terms (y ^ v 1 ). The result in free space is recovered by taking the limit ro << T\ (i.e. very distant wall). Similar 
consideration holds for the other contributions we are now going to evaluate. 



In the continuum limit, the contribution from second-order fields in (|T5[) becomes 



(0 kj ,U |E(°)(r,i) -E( 2 )(r,t)|0 kj ,U) =-^E E - r„) 



/ dk 

Jo 



— ikct 



VP 



k + kn 



pqq' vv'— 0,1 



, — ik {ct—r v ) _ ik(ct—r u )\ 



+ 



kn — k 



o ife (ct— Tv) _ ik(ct-r v )\ 



F, 



sin(fc7v) 



(20) 



7 



Finally, after some straightforward algebra, we obtain 

(Okj,U |E(°)(r,t).E( 2 )(r,f)|0 kj -,U) + (0 k ,,U |E® (r, t) ■ E<°» (r, t)\0 kj , | A ) 

1 1 r+oo 

x [■^^(fco 3 /p 9 (fco^)/^(fcr I ,Oe- 4fco(ct -^ ) e^ fc ( ct ^') + 

+ feg/ M (*br„)/ M /(-*r^)e- < * ( ct - r ")e- tt t ct+r -'5 + fc 3 / P9 (-^)/ P9 '(-^0e~ Jfc(r " +r -' ) ) 
+ ^ I (-fc 3 / P9 (-fcor,)/ P9 '(fcr„Oe 4 ' £o(ct -^ ) e- 4fc ( ct -^') 

- ^/ M (-feor w )/^(-fc^')e ifeo(ct - r " ) e- ife(ci+r >'' ) + k 3 f pq (-kr„)f pq , (-kr v ,)e- ik ^ +r ^ 



(21) 



c.c. 



where the terms with a pole in k — fco should be taken 
as their principal value. We now extend the integral over 
k to the complex plane, closing the integration path at 
infinity in the fourth quadrant and on the negative imagi- 
nary axis. Assuming no other singularity in the integrand 
function except the pole at fc = fco and that it vanishes 
at infinity, the result can be written as an integral along 
the imaginary axis plus the pole contribution 

poo poo 

P I dkg{k) = —i I dk g(—ik) — iirRes{g(k)}k = k a 
Jo Jo 

where g(k) indicates the integrand function. 



We now give the explicit expression of the expectation 
value of the square of the electric field for t > r± /c, that is 
after that the direct and reflected waves have reached the 
observation point (this quantity vanishes for t < ro / c, of 
course, but it is not zero for r$/c < t < r±/c due to the 
"direct" contribution). In this timescale, the contribu- 
tion from the pole at fc = fco for t > T\jc in (f2"2")) exactly 
cancels the contribution (fill)]) when the three second-order 
contributions in (|15p are summed up, exactly as in the 
case of an atom in the free space discussed in [l2j]. We 
obtain 



(o kj ,u is 2 (r,t)io kJ ,u)-< = --E E {-lro-'-H-ifo-'" 

pqq 1/1/—O 



dkk 6 a qq >(ik)f pq (ikr„)f pq ,(ikr„,)e fe ( r -+ r v) 



^° / dkk 3 -- kct 



(zfc)(fc Re f pq (k r„)e-^ ct - r ^ - fc Im / p? (fc r,,) e - ifc °( ct ~^) 
(f pq '(-ikr u ,)e kr »' + f pq '{ikr v >)e- hr "'\ 



(23) 



r 



where we have explicitly indicated that the "zero-point" 
electric energy density E^ p due to zero-order field oper- 
ators has been subtracted, and we have also introduced 
the dynamic atomic polarizability tensor at imaginary 
frequencies 



a qq ,(ik) = 



2fc, 



hc(k 2 



(24) 



It is worth considering the limiting case of Eq. 
for t — > +oo. In this limit, fast oscillating terms in- 
side the fc-integral average to zero, finally yielding the 
time-independent value given by the first two lines only 
of Eq. (|23|) . For an isotropic atom or polarizable body 
(Px = (J* y = /i z ) this becomes 
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(0 kj ,lA\E 2 (r)\0 kl ,U) - Eg = —J2 / *fc 6 a(*fc)(-rrr 

n ^ v 
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7T 



fcr 
1 



^3^3 ^4^.4 £-6^6 

+ (<5pg - 3f 0p f 0(J 



-2fcr„ 



(25) 



-fc(r +ri) 



where a(fc) is the isotropic electric polarizability of the 
atom. There are three different terms that contribute 
to the electric energy density. The first comes from the 
v = v' = contribution and it is identical to that found 
in for a two-level atom in the free space. The second 
term is the v = v' = 1 contributions and has the same 
functional expression of the previous one with ro replaced 
by t*i (distance of the observation point from the image 
of the atom). Finally, there is a mixed term depending 
on both distances ro and r\. The physical meaning and 
role on interatomic Casimir-Polder interactions of these 



terms will be discussed in detail in the next Section. 

Similarly to the electric case, from the results of Sec- 
tion IIIII we can calculate the magnetic energy density, 
which we shall use in the next Section in order to calcu- 
late magnetic interatomic Casimir forces in the far zone. 
We shall not give the explicit expressions of the time- 
dependent magnetic energy density, but only its asymp- 
totic expression for t — ► +oo. In this limit it settles to 
a time-independent value, as the electric energy density, 
given by 



(0 kj ,U|H 2 (r)|0 k „U)-< = ^ £ £ 

r+oo £,6 

x / dk 2 

JO Kq + K 



p,q,q' v^v'— 0,1 



dk TTT^9 Pq {ikr l ,)g pq ,(ikr v ,)e- k{r » +r ^ (26) 



p.q.q' v.v' — 0, 1 



dk k^a qq \ik)g vq (ikr v )g m > {ikr v >)e~ k{ - ru+r » ,) 



(similarly to the electric case, the zero-point magnetic energy density has been subtracted) where 



(ifcr)=^e Ms f s (^ + p LY (27 ) 



For an isotropic molecule it becomes 

/•+OC 

dkk 6 a{ik)g pq {ikr u )g pq {ikr u ,)e- k{r » +r »' ) . (28) 



(0 k „U|i? 2 (r)|0 kj ,U)-< = ^ £ / 

p,q v,v<=0,l 



V. THE FAR-ZONE ELECTRIC AND 
MAGNETIC INTERATOMIC CASIMIR-POLDER 
ENERGIES WITH BOUNDARY CONDITIONS 

In this Section we apply the results of the previous 
Sections for the time-dependent Maxwell fields to the 
calculation of dynamical electric and magnetic Casimir- 
Polder interactions between two atoms in the presence of 
the conducting wall, also recovering known results in the 
asymptotic regime. The electric and magnetic field oper- 



ators of Section Hm and the related expectation values of 
the energy densities obtained in Section HVl can be used 
to obtain static and dynamical Casimir-Polder forces be- 
tween two atoms in the presence of a boundary condi- 
tion such as a conducting wall. In fact, the interatomic 
Casimir-Polder interaction in the so-called far zone, that 
is for distances between the atoms larger than their typ- 
ical transition wavelength from the ground state, can be 
obtained from the interaction of one atom, described as 
a polarizable body, with the field fluctuations generated 
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by the other atom [2l], [22|, [23| ■ This approach has been 
recently proved to be useful also in the case of dynamical 
(time-dependent) Casimir-Polder forces for atoms in ex- 
cited states [13] or in partially dressed states [l3j]. Thus 
we evaluate the (electric) Casimir-Polder interaction en- 
ergy in the far zone between our "source" atom A and a 
second atom B placed in r E , when the conducting wall 
in present, as 



AE F 



■-a B E (E 2 {v Bl t)) A 



(29) 



where a E is the static electric polarizability of atom B, 
assumed isotropic for simplicity, and the average value 
of E 2 (r,t), evaluated at the position of atom B, is cal- 
culated on the initial state of the field and of atom A 
(for example the bare ground state, as in the previous 
Section). In a quasi-static approach, the Casimir-Polder 
force is then obtained as minus the derivative of the in- 
teraction energy with respect to the distance between 
the two atoms. Effects related to the atomic motion, as 
discussed in [H,[I3i are neglected here. Besides being 
well suited for dealing with dynamical interactions, our 
method also allows a simpler derivation of static Casimir- 
Polder forces and a transparent physical interpretation, 
compared to the usual methods based on fourth-order 
perturbation calculations [24| . 

From our results dynamical magnetic Casimir-Polder 
interactions can also be easily obtained, such as the inter- 
action of a magnetically polarizable body B (with static 
magnetic polarizability a E ) with the magnetic field fluc- 
tuations generated by atom A, as 



AE M = -~< 



[H 2 {v B .t)) J 



(30) 



The relations above show that the Casimir-Polder in- 
teraction energy between two atoms is proportional to 
the expectation value of the square of the fields generated 
by one of the two atoms at the position of the other atom, 
that is to the average electric (and magnetic) energy den- 
sity or equivalently to the field fluctuations generated by 
the atom. 

In order to evaluate equation (|29[) in the far zone, we 
can use the expression (1231) of the electric energy density 
with r = vb (in the limit t — > oo), approximated by sub- 
stituting the dynamical polarizability of the atom with 
the static one; this allows to perform the k integral. We 
finally obtain 



1 



23hc 



AE E = -~ a «{E i {v B )) A = - ^ 



y E a E 



1 



1 



+ 



16hca E a E 



7iTo?i(ro + ri) E 



3r r x + p (r + 5r ri + r x ) 



with 



p = 7'o sin Oo — r\ sin f/i 



(31) 



(32) 



where 9o (#i) is the angle between the vector ro (ri) and 
the perpendicular to the wall. 

Expression (|3ip . which was already obtained with dif- 
ferent methods in [3J, [J], shows that in the presence of 
the conducting plate, the Casimir-Polder interaction be- 
tween the two atoms A and B consists of three terms: the 
r~Q 7 interaction between the two atoms as in the absence 
of the wall; the r^ 7 interaction between one atom and 
the image of the other atom as reflected by the wall; a 
term involving both coordinates r and r\ . The first two 
terms give an attractive component, whereas the third 
onegives a repulsive component. However, as discussed 
in [4| , the total interaction is attractive for any geometric 
configuration of the atoms with respect to the wall. 

Our results in the previous Section also allow to eval- 
uate the time-dependent Casimir-Polder interaction be- 
tween the two atoms, using the time-dependent expres- 
sion for the squared electric energy density (|23p (and not 
only its asymptotic value for t — > oo). The interaction en- 
ergy between the two atoms in the far zone as a function 
of time for ct > r x , that is after that both the direct and 
reflected fields generated by atom A have reached atom 
B, is characterized by oscillations around its asymptotic 
value, with an amplitude decaying with time to zero. For 
ct 3> r\ the energy settles to its stationary value, as ob- 
tained in a time- independent approach [j| or from (|31jl . 
This behavior is analogous to that already found in the 
time-dependent Casimir-Polder force between an atom 
and a conducting wall [l6j |. 

In a similar way we can explicitly evaluate the mag- 
netic far-zone Casimir-Polder interaction energy, in the 
presence of the wall, by substituting (|28|) into (|30| . In the 
far zone, after substitution of the dynamical polarizabil- 
ity of the atom with the static one and the k integration 
we get 



AE M = -\ a 



i ir 7 he 

,H {r B ))A = -^—oc^a 



A B 
E a M 



1 



1 



16hca E af 4 ro cos #o r i cos 6\ 



■Kr%r\{ro + nY 



"o + 5r o^i + r i ■ 
(33) 

We note that, similarly to the electric case, the mag- 
netic interaction contains three terms: two are related 
to the atom-atom (depending from ro only) and atom- 
image (depending from r\ only) interactions, and the 
other one depends on both coordinates ro and T\. It 
should be also noted that the magnetic Casimir-Polder 
interaction between the two atoms is repulsive. How- 
ever, in ordinary situations magnetic Casimir-Polder in- 
teractions are much weaker than electric ones, the to- 
tal force being attractive. For large atom-wall distance 
[r\ ^> ro), Eq. (f3"3"|) reduces to the well-know station- 
ary electric-magnetic Casimir-Polder interaction between 
two atoms in free space [111 [26|. Finally, we wish to 
stress that the atom-atom dynamical Casimir-Polder in- 
teractions discussed in this Section add, of course, to 
the know n dy namical atom- wall interaction energies, dis- 
cussed in [lj| [l6|, [ItJ for the case of electric interactions. 
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VI. CONCLUSIONS 

In this paper we have calculated the time-dependent 
electric and magnetic Maxwell field operators in the 
Heisenberg representation for a two-level atomic source 
near a boundary condition such as a perfectly conducting 
wall. The method used consists of an iterative solution 
of the Heisenberg equations for the field operators. The 
results obtained for the field operators are quite general, 
not depending on the specific atom-field initial state. We 
have explicitly obtained solutions up to the second order 
in the electric charge and showed that they can be sep- 
arated into a direct and reflected term. The latter is 
related to the presence of the wall. Both terms propa- 
gate in a strictly causal way. We have then evaluated 
the time-dependent electric and magnetic energy densi- 
ties around the atom for a specific initial state of the 
system, that is the atom in its bare ground state and the 
field in its vacuum state. This has finally allowed us to 
evaluate the dynamical (time-dependent) Casimir-Poldcr 
interaction energy between the "source" atom and a sec- 
ond atom, considered as a polarizable body, when both 
are placed in the vicinity of the conducting wall. Both 
electric and magnetic Casimir-Polder energies have been 
considered. We have also shown that for t — > oo both 
the electric and magnetic dynamical Casimir-Polder en- 
ergies settle to time-independent values, which coincide 
with those that can be obtained by a time-independent 
approach. We stress that our method is general and not 
limited to a specific initial state of the system; in fact 
our solutions for the field operators in the presence of 
the wall can be used for any atom-field initial state, be- 
ing operator solutions. 



APPENDIX: POLARIZATION SUM 



In this Appendix we evaluate the following expressions 
which have been used in Section Ell 



£ f dO(f ki (r)) p (f k ;(r A )), 
j J 

£ J dfi(V x f kj -(r)) p (f kj -(r A )) g 



(A.l) 



where p,q — x, y, z. In order to show the method used, 
we first make a specific choice of p and q, namely p = x 
and q = y, and we calculate the first expression of (|A.1|> . 
The technique is the same for all the other cases. Using 
the mode functions ((2]) and applying the sum rule for the 
polarization unit vectors 



k h 



(A.2) 



where 6 pq is the Kronecker delta function, we get 
J dQ(f kj (r)) x (f kj (r A )) y 

cos(k y L/2) sin(fc z z) sin(fc z G?). 



■ d^~i krpk lt cos 

k 



x sin 



(A.3) 



We can now exploit the fact that in the limit L — > +oo, 
the oscillating functions containing L can be replaced by 
their mean value over a period; for example 



sin 2 (fc x Z/2) -> 1/2 



sm(k x L/2) cos(fe x £/2) -> 0. 



(A.4) 



Thus we get 



k 2 
2 

k~ 2 



sin(k x x) sm(k y y) sin(fc z z) sin(fc z d) 

cos(k x x) cos(k y y) sin(fc z z) cos(fc z rf) 



(A.5) 



Finally, using 



sin(kD) 

dfl cos(k x a) cos(kyb) cos(k z c) = 4n — (A. 6) 



with D — y/a 2 + b 2 + c 2 , we get 



Y, [ dO(f kj (r)) x (f kj (r A )) H = ^Y1 V * V ? 



sin(fcr„) 



(A.7) 



where ro = y 1 x 2 + y 2 + (z — d) 2 and r\ = 
\J x 2 + y 2 + (z + d) 2 . All the other terms can be 
obtained using the same method. The general results 
are 



I dfi(f ki (r)yf ki M), 

3 

- 4ir V/_ii"(i^)p sin ( fcr ^) 



(A. 



v=0 



and 



I *(Vxf kj (r)) p (f kj (r,)), 

3 J 



(A.9) 



r/=0 



where F pq and G pq are the differential operators defined 



m Section [ml 
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